Sea quark effects in B K from Nf = 2 clover-improved Wilson fermions 
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We report calculations of B K using two flavours of dynamical clover-improved Wilson lattice fermions and look for depen- 
dence on the dynamical quark mass at fixed lattice spacing. We see some evidence for dynamical quark effects. In particular 
B K decreases as the sea quark masses are reduced towards the up/down quark mass. Our meson masses are quite heavy and a 
firm prediction of B K is a task for future simulations. 



1. INTRODUCTION 

Bk determines the value of the AS = 2 neutral kaon 
mixing matrix element according to 

( fO|gAS=2 0u)| ^ : O ) = % BK(M)f 2 m 2^ (1) 

with fi indicating the scale dependence of the operator 
e AS = 2 0u) = ay„(l - y 5 )d - y 5 )d. 

The quenched lattice value of Bk has more or less 
settled down (see (T] for a recent review) and the 
1997 staggered fermion result fi^(MS, 2 GeV) = 
0.63(4), corresponding to a renormalisation group in- 
variant (RGI) value Bk - 0.87(6), remains the bench- 
mark. The quoted error does not include quenching 
effects which have been estimated to be as high as 
+ 15% |3| and remain the primary systematic effect 
to be addressed. Here we report on a calculation |4| 
using two degenerate flavours of dynamical (clover- 
improved) Wilson fermions, aiming to look for sea- 
quark dependence in Bk- 

2. SETUP AND SIMULATION 

The operator of interest in the continuum is 
QiQi) = Ijnd 1-fd + Ijujid Jy^ysd, (2) 
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which is the parity conserving part of Q AS=2 (fi) in 
eq. Q. For Wilson fermions, owing to the explicit 
breaking of chiral symmetry, this operator mixes with 
other four-fermion operators. Including the overall 
multiplicative renormalisation, the subtraction of the 
unwanted operators may be expressed as 

QiQi) = ZQx, g\) (<2' att + Yj A/(<7o)Gl att )- (3) 

The coefficients Z and A, have been determined per- 
turbatively to 0(a s ) for MS-NDR in 06|. The AS = 
2 lattice operator can now be matched to the con- 
tinuum one, having the same chiral properties up to 
O(a ). There can be remnant operator mixing, begin- 
ning at 0(a 2 a°), because the A,- coefficients are not 
known exactly. Even though we use an improved- 
clover action, the four-fermion operators are unim- 
proved and O(a) artefacts can also be present. 

By calculating matrix elements with non-zero mo- 
mentum kaons we introduce another degree of free- 
dom, allowing us partially to remove the lattice arte- 
facts. The chiral behaviour of the matrix element can 
be written as 

</ i0 ,p|£i(A')|.P ( \q> = a' +B'm 2 P + S'm 4 P + 

p-q [y + y' + (e + e')m 2 P + (f + ?)p-q\ + ■■■ (4) 

The parameters y', e' and are corrections to the cor- 
responding physical contributions. In contrast, a', B' 
and 8' are artefacts of chiral symmetry breaking, in 
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Table 1 

The configurations used |8 9|. Lattice spacings are 
fixed from ro. mp/my is quoted for K sea = K m \- 



p 


csw 


'''sea 


a/fm 


mp/my 


#cfgs 


5.20 


2.0171 


0.1350 


0.103(2) 


0.70(1) 


100 


5.26 


1.9497 


0.1345 


0.104(1) 


0.78(1) 


100 


5.29 


1.9192 


0.1340 


0.102(2) 


0.83(1) 


80 



this case 0(a, a 2 a°) with a 2 small, which have to be 
subtracted. In particular the a' term makes B K di- 
verge in the chiral limit. Following |3|, we ignore 
chiral logarithmic terms in eq. (0}, given the errors 
and large masses used here. For our calculation we 
neglect higher order terms and use the following ex- 
pression for the matrix elements: 

(P°,pieiOu)|f ,q) = a' + B'm 2 p + {j + y')p-q. (5) 

We calculate Bk using Clover-improved Wilson 
fermions 1 7 1 on the UKQCD set of unquenched con- 
figurations 1 8 9 1 listed in tabled We use three differ- 
ent sea quark masses in the region nip/my > 0.7 on a 
volume of 16 3 x 32 (mpL > 7) with a nearly constant 
lattice spacing determined from the Sommer scale, ro. 
Propagators and correlators are calculated using the 
FermiQCD 1 1 11 1 code. Five valence quark prop- 
agators at k = 0.1356, 0.1350, 0.1345, 0.1340 and 
0.1335 are generated for each sea quark. We follow 
the standard procedure 1 12 1 to evaluate the matrix el- 
ement from 3- and 2-point correlation functions. 

C (y) {t x ,t y ;p x , PlJ ;p) = 

J^sfc f -)G(°; t y )}e^ e*» y , (6) 



(7) 



We use local pseudoscalar sources 7, for the kaons 
(smearing with the method of 1131 did not improve 
the signal significantly). In the 3-point functions, Q is 
at the origin and t y is fixed at 10 while t x is varied over 
the full temporal range. We have checked neighbour- 
ing values of t y but observe no dependence, imply- 
ing that the ground state is reasonably well isolated. 
The momentum configurations {p A , p^} are chosen 
from the pairs {(0, 0, 0), (0, 0, 0)}, {(0, 0,0), (1, 0, 0)} 
and {(1, 0, 0), (0, 0, 0)} where an average over equiv- 
alent configurations is implied. The matching co- 



efficients are evaluated using the boosted coupling 
gi = 6(P)/B, where the average plaquette values are 
(P) e {0.5336, 0.5399, 0.5424} on our ensembles. 



3. ANALYSIS AND DISCUSSION 

The matrix element is extracted using the ratios: 
C {3) (t x ,t u ; p x ,p„;u) 

* 3 (p) = 2 (2) ; 2) (8) 



X(0) = - 



C {2) (t ) 



Cf p {t x ) 



X(jt) = X(0) 



PxPy 



(9) 



where Za is the axial current renormalisation. Our 
Method I is to follow 1141151 and fit the equation 



R 3 (p) = a' +B'X(0) + (y + y' )X(p), 



(10) 



to find B K from y, neglecting y' . This gives an es- 
timate of the leading term in an expansion of B K for 
each set of different valence quarks with a given sea 
quark mass, with the kaons not necessarily being at 
the physical kaon mass. Using a linear fit versus the 
unitary pseudoscalar masses, (amp) 2 \ Km=K ^, to go to 
the up/down limit gives 

B K (MS, 2 GeV) = 0.49(13), (11) 

corresponding to the RGI value Bk = 0.69(18). In 
this method the pseudoscalar masses are for valence 
quarks in the simulated region {m q w 2m s ) and the 
B K estimate is one where the sea quarks are realisti- 
cally light but the valence quarks are heavier than the 
physical strange quark. 

In Method II we follow |5 16 1 and estimate Bk for 
each (/e S ea, A'vai) combination: 



/? 3 (p)-fl 3 (0) 



V(A = B K (p,K se:i ,K ml ). (12) 

x(p) - x(0) (lWvui) 

The resulting values are illustrated in fig. E Taking 
the unitary points (as done in 1171 ~). with /f sea = A- va i, 
and extrapolating to the physical kaon mass results in, 



B/f(MS,2GeV) = 0.48(13). 



(13) 



m^ hys , but the sea and valence 



Here we have m P 
quarks are degenerate and hence the sea content is not 
as light as the up/down quarks. Moreover, the kaon 
comprises two quarks around m s /2. 
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Figure 1. Z?A-(MS,2GeV) for each (A" se a, Kvai) as a 
function of the squared pseudoscalar mass. Dashed 
lines are for guidance to distinguish sets with differ- 
ent sea quarks. The filled points joined by a solid line 
have A: sea = K m \. The lightest point (marked by a large 
cross) is excluded from the analysis. 



Although we recognise the presence of lattice arte- 
facts, it seems that dynamical quark effects can lower 
the value of Bk- Taking this with the observation 1 18 1 
that Nf =2 results are always below those for Nj — 0, 
a statement also valid for subsequent works, we see 
that when one has two still-heavy sea quarks, Bk 
starts to decrease but is consistent with the quenched 
value. When the sea quarks are taken towards the 
massless limit, Bk becomes distinctly lower. A simi- 
lar effect is seen in the domain-wall results of 1 1 71191 . 
In contrast, preliminary 2+1 flavour results with im- 
proved staggered quarks 1 20 1 show some sign of a de- 
crease in Bk with decreasing sea-quark mass, but the 
overall result is consistent with the quenched value. 
It is intriguing to note that recent Unitarity Trian- 
gle fits, which determine Bk after imposing the re- 
maining constraints, also find a value below the usual 
quenched lattice result (for example, B K = 0.65(10) 
in 1121 1 ). However, our primary result is a decrease in 
Bk as the sea-quark mass decreases, rather than a firm 
value for use in phenomenology. 



5. 



Acknowledgment: we thank Ivan Wolton, Ian Hardy, 
Oz Parchment and the Iridis team at the University of 
Southampton for computing support. 

REFERENCES 



R. Gupta (2003), hep-lat/0303010. Contributed 

to Workshop on the CKM Unitarity Triangle, 

CERN, Geneva, Feb 2002. 

JLQCD Collaboration, S. Aoki et al., Phys. Rev. 

Lett. 80 (1998) 5271, hep-lat/9710073. 

S.R. Sharpe, Nucl. Phys. Proc. Suppl. 53 (1997) 

181,hep-lat/9609029. 

UKQCD Collaboration, J.M. Flynn, F. Mescia 

and A.S.B. Tariq (2004), hep-lat/0406013. 

R. Gupta, T. Bhattacharya and S.R. Sharpe, Phys. 

Rev. D55 (1997) 4036, hep-lat/961 1023. 

S. Capitani et al., Nucl. Phys. Proc. Suppl. 63 

(1998) 874, hep-lat/9709049. 

M. Luscher et al., Nucl. Phys. B491 (1997) 323, 

hep-lat/9609035. 

UKQCD Collaboration, A.C. Irving, Nucl. Phys. 

Proc. Suppl. 94 (2001) 242, hep-lat/0010012. 

UKQCD Collaboration, C.R. Allton et al., Phys. 

Rev. D65 (2002) 054502, hep-lat/0 107021. 

M. Di Pien-o, Nucl. Phys. Proc. Suppl. 106 (2002) 

1034,hep-lat/0110116. 

M. Di Pierro (2000), hep-lat/001 1083. 

M.B. Gavela et al., Nucl. Phys. B306 (1988) 677. 

UKQCD Collaboration, P. Boyle, J. Comput. 

Phys. 179 (2002) 349, hep-lat/9903033. 

SPQCDR Collaboration, D. Becirevic et al., 

Nucl. Phys. Proc. Suppl. 119 (2003) 359, hep- 

lat/0209136. 

M. Crisafulli et al., Phys. Lett. B369 (1996) 325, 
hep-lat/9509029. 

R. Gupta et al., Phys. Rev. D47 (1993) 51 13, hep- 
lat/9210018. 

T. Izubuchi (2003), hep- 



10. 

11. 
12. 
13. 

14. 



15. 



16. 



17. 



18. 



19. 



20. 



21. 



RBC Collaboration, 
lat/03 10058. 

A. Soni, Nucl. Phys. Proc. Suppl. 47 (1996) 43, 
hep-lat/95 10036. 

RBC Collaboration, C. Dawson (2004). These 
proceedings. 

HPQCD-UKQCD Collaboration, E. Gamiz et al. 
(2004), hep-lat/0409049. These proceedings. 
UTfit Collaboration, M. Bona et al. (2004), hep- 
ph/0408079. 



